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Abstract 



lO ■ Within tlie framework of Bolimian meclianics dwell times find a straight- 

xl ' forward formulation. The computation of associated probabilities and 

distributions however needs the explicit knowledge of a relevant sample of 
trajectories and therefore implies formidable numerical effort. Here a tra- 
jectory free formulation for the average transmission and reflection dwell 



o 
\o 
o 

j^— s . times within static spatial intervals [a, b] is given for one-dimensional scat- 

■""•^ ' tering problems. This formulation reduces the computation time to less 

i-Q ' than 5% of the computation time by means of trajectory sampling. 

^ ! PACS numbers: 03.65.-w 

O-^! 1 Introduction 



For a ID static detector located in the spatial interval [a,b] (see figure 1), the 
average dwell time of an ensemble of quantum systems vifith wave function ^ 
from time tj up to time tj is generally agreed to be given by 

^' dt f d^ |vl/(i,0l', (1) 

i ^ a 

cf the reviews [1]^. It is motivated by classical reasoning [2], and it also has been 
derived within Bohmian mechanics [3]. In a recent work [4] a corresponding 
dwell time operator has been investigated. 

Dwell times of the type (1) were associated with interaction times in collision 
processes long ago [5] . The relevance of these interaction times to time resolved 
scattering experiments has been studied, e.g., in [6]. Through these works it be- 
came clear, that differences of average dwell times formed between a scattering 
wave packet and its free incoming asymptote are measurable quantities. 

1 Clearly expression (1) does not exist for bound states. 



With the introduction of the Larmor clock by Baz [7] , dwell time expressions 
themselves have got measurable status. The idea is that a small and uniform 
magnetic field, which is confined to a small region of space, causes a Larmor- 
precession of the spin-polarization vector of the scattered wave. It was shown 
in [8] that the Larmor clock indeed reveals the average dwell time. If it also is 
capable of displaying the spectral distribution of some dwell time operator is 
still unclear. 

The specialization of the Larmor clock to the case of one-dimensional scat- 
tering was done by Rybachenko in [9]. In this work (to the authors knowledge 
for the first time) a distinction into the dwell times of the finally transmit- 
ted respectively reflected partial waves has been introduced. For short such 
selective dwell times will further on be denoted as transmission and reflection 
times. Another approach to transmission and reflection times, grounded on a 
specific experimental scheme, was introduced by the oscillating barrier model of 
Biittiker and Landauer [10]. This model is widely believed to have ignited the 
tunnelling time controversy anew. A further extension of the Larmor clock by 
Biittiker [11] incorporates the effect of spin-alignment with the magnetic field. 

More recently the interest in transmission times has been driven on the 
one hand by the very indirect time measurement techniques of the condensed 
matter community, especially in connection with tunnelling in semiconductor 
heterostructures or Joscphson junctions. On the other hand, the progress in 
laser cooling techniques in quantum optics delivers another valuable tool for 
future time resolved scattering experiments. 

It is no surprise that the detailed definition of transmission and reflection 
times depends on the situation under scrutiny. A systematic operator approach 
embodying several such possible deflnitions was given by Brouard et al [12]. 
It was shown, however, that transmission and reflection times derived within 
the framework of Bohmian mechanics arc not included in this catalogue ([13], 
chapter 5). 

Bohmian mechanics comprises the mathematical concept of world lines or 
trajectories and with this the term 'particle' obtains substance in quantum the- 
ory again. Thus the notion of dwell time can be addressed in a straightforward 
manner, very much like in classical mechanics. The dwell time of a particle 
in the spatial interval [a, b] is simply defined as the duration during which the 
particle's trajectory is localized within [a, 6]. For an elaborate discussion of 
Bohmian mechanics see [14]. 

As the numerical effort involved with the calculation of Bohmian world lines 
is immense, there have been attempts to compute the Bohmian transmission and 
reflection times without the sampling of trajectories. A related one-dimensional 
bound-state situation has, e.g., been studied by Stomphorst in [15]. In the 
present paper a formulation without trajectories for the transmission and re- 
flection times in a genuine scattering situation is presented. The derivation 
closely follows ideas developed for the treatment of ID arrival time in [16]. As 
an application the scattering from a double potential barrier, reminiscent of 
semiconductor-hetcrojunctions in, e.g., resonant tunnelling diodes, is consid- 
ered. 



2 Bohmian transmission and reflection times in 
ID scattering situations 

Wc consider a scattering situation in which the M0ller operators fi™ and ri°"* 
exist and are asymptoticaUy complete [17]. *(i, •) = e"'^'/'^l^™0o = Q'"'e~'""*^'^(, 
denotes the scattering solution with incoming asymptote (t){t,-) ~ e^'^"*'''(/)o. 
The solution <j>{t, •) of the free Schrodinger equation is chosen to be localized on 
the negative spatial semi-axis for < ^ — oo. That is the case if and only if the 
Fourier transform J-'{<f>o) is localized on the positive half-line [18]. 

The Bohmian transmission time in the above scattering situation is defined 
as follows. Let ^x{t) denote the Bohmian trajectory which at time t = passes 
through the point x. There exists a critical value Xc G K, such that lim jxit) = 

t — >oo 

— oo for all X < Xc and that lim 7x(i) = oo for all x > Xc- Therefore the 

t — »oo 

Bohmian transmission time is defined as 

{tt)= J dx\^{0,x)\''-eix~Xc)Jdtx[aMhAt)) (2) 

— oo Ti 

with X[a,b] the characteristic function on the interval [a, b] and © the Heaviside 
step function. See, e.g., [19]. Analogously the reflection time < r/f > is defined 
by replacing the term Q{x — Xc) in the right hand side of equation (2) by 
6(.Tc — x). The critical trajectory ^x^it) (71;^ (0) = Xc) is implicitly defined by 

/•oo 

|r|2= / \^{t,0\'d^, VteM. (3) 

"'7x,(t) 

Thereby 

/>oo 

|T|2:= lim / |*(t,OI'dC (4) 

t— oo_/q 

is the transmission probability of the scattering system. The lower limit of the 
integral in (4) can equally be replaced by any finite g £ K. Accordingly by 
|i?|^ := 1 — |T|^, the reflection probability is defined. Thus the conditional 
transmission respectively reflection times, i.e. transmission and reflection times 
normalized to the fraction of transmitted respectively reflected particles of the 
entire ensemble, arc < tt >'^:— jTinT < tt > and < tji >'^:= yirj < th >. 



From the computational viewpoint, the two terms < tt > and < r/j > by means 
of Bohmian mechanics are achieved in a straightforward manner. One chooses 
an appropriate sample of initial values on the configuration space Eg = M, cal- 
culates the corresponding trajectories over a sufficient range of time, determines 
the dwell time for each trajectory, labels the trajectories as transmitted or re- 
flected according to their position at large times and finally calculates, according 
to the weight of each trajectory and just as in classical statistics, the average 
times. That this programm involves formidable numerical effort is evident. 



G 
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Figure 1: Spacetime region G = [a,b] X [Ti,Tf]. 

However, the calculation of Bohmian transmission and reflection times can 
be reduced to the computation of current density integrals along the edges at 
X = a and x = b. The next proposition represents a generalization of expressions 
already proposed in [20]. 

Proposition 1. For ID scattering solutions ^^(pt with Q{K)(J)q = (I>q, \\(I>o\\ = 
1, for the Bohmian transmission and reflection times within [a,b] x [r^jTy], hold 

{tt) - r dt [min {fait), |Tp} - min {f,{t), |Tp}] (5) 

and 

(tr) = r dt [max{/,(i),|T|2} -max{/b(f),|Tp}] (6) 

with 

fq{t) := / Jis,q) ds. 

Here j is the quantum mechanical probability current density. 

An essential ingredient for the proof of proposition 1 is the relation 

/t /'CK3 

j{s,q)ds= I'^it^Ol'dt (7) 

-oo J q 

It depicts that the probability of finding a particle to the right of q at time t 
is equal to the amount of probability, which has passed q up to the time t. A 
plausibility argument for equation (7) and a rigorous proof for a limited class of 
scattering situations is given in appendix A. The proof of proposition 1 is given 
in appendix B. 

Remark: The formulation of Oriels et al [20] assumes the case in which the 
current density at the right edge of the barrier doesn 't change its sign and is 
positive for all times. In this case fb{t) < |Tp, Vi G M, because 



IT|2 (4) r 
\1 \ ~ lim 

t — 'OO 



oo t oo 

J |vl/(i,0r d^ ^^ ^Inn j j{s,b) ds= J j{s,b) ds 



and further 



Mt)^ I j{s,b)ds'^''T'' I j{s,b)ds^\T\ 



Then equation (5) reduces to 



■^f 



dt 



\Ut),\T\^] - Ut)] 



and equation (6) to 

^^"^^ ^W^^'' ^^Wlr '^^ [max{/„(t),|rp} - |T|2] 

which reproduces equations (14) and (15) of [20] for the special choices n = 

and Tf = oo. 

Obviously an interesting task would be to construct a device, i.e. a clock, which 
measures the Bohmian transmission and reflection times. Such a clock should 
display the respective time, let us say, through the Bohmian center of mass posi- 
tion of its hand at the instant of its readout, which presumably has to be chosen 
by the experimenter. The combined system's wave function would be modelled 
by an appropriate Schrodinger equation, incorporating the interaction between 
the micro-system and the clock. There is a self-adjoint operator corresponding 
to the hand's positions. The crucial question is whether this observable's spec- 
tral distribution in a given state at the instant of its readout coincides with the 
respective Bohmian transmission (reflection) time distribution. Probably this is 
not the case for all states. However, similarly to the issue of exit time statistics 
[21], a subspace might be identified on which the Bohmian distribution coincides 
with one of the standard quantum mechanical distributions. 



3 Numerical example: transmission and reflec- 
tion at a double potential barrier structure 

As an example for a ID scattering situation we consider the case of a Gaussian 
wave packet impinging on a double potential barrier, i.e. we are looking for 
solutions '^ to the Schrodinger equation 



ihdt'S 



2m 



dl + V 



* 



with 



V = Vo- X[a,b] + Vl ■ {X[a',a[ + X]b,b']) ■ 

X[a,f3] denotes the characteristic function on the interval [a, /3] G M and a' < a < 
b < b' (see figure 2). The parameter reduction h,m,Vo — » 1, e.g., is achieved 



by taking time in units of 2^, space in units of 



Vi is taken in units 



of Vq. In figure 3(a) the evolution of the probabihty density |^p is given for 



the case 



-6, a 



-3, 6 = 3, 5' = 6 and V^i = 2 in the chosen units. The 



mean kinetic energy of the packet is l.S^Vb = 2.25Vb. In figure 3(b) a sample of 
50 corresponding Bohmian trajectories is illustrated. The initial distribution of 
starting points of the trajectories resemble the initial Gaussian distribution of 
the wave packet. Figure 4 shows a zoom into the area indicated by the rectangle 
in figure 3(b). 



a' a b b' ) 

Figure 2: Double potential barrier. 
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Figure 3: (a) Evolution of the probability density I'J'p. (b) Corresponding 
sample of 50 Bohmian trajectories. The double potential barrier 
is indicated by the hatched area. 



In figure 4 it becomes clear that, as trajectories change their direction at x = 6, 
the current density in this case changes its sign also at the right edge of the area 
in question. Therefore the restricted formulae of Oriols et al loose their validity 
and the generalized expressions (5) and (6) have to be applied. 
Figure 5 shows the current densities and corresponding integrated current den- 




Figure 4: Zoom into the region indicated in figure 3(b). 







Figure 5: Current densities j(-,a) (a) and j{-,b) (b) and corresponding in- 
tegrated current densities /a(-) (c) and /(,(■) (d). 



sities at the edges x = a and x = b, respectively, as a function of time. In (c) 
and (d), the transmission coefficient jTp = 38.31% is indicated by the dashed 
Hne. 

Finally, in figure 6 the mappings Tt{s) (a), Tr{s) (b) and Td{s) (c) are 
illustrated, which give the Bohmian transmission, reflection and overall aver- 
age dwell times < tx >, X G {T,R,D} in [a, 6] from time Ti — Q onwards 
as a function of the upper temporal bound s — Tf. In addition the conditional 
transmission and reflection times, i.e. normalized to the fraction of finally trans- 
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Figure 6: (a) Transmission times TrC*) (solid line) and corresponding con- 
ditional transmission times (dashed line), (b) Reflection times 
Tfi(t) (solid line) and conditional reflection times (dashed line). 
(c) Average dwell times Tb(t) inside [a,b] X [0, t]. 



mitted or reflected particles, are indicated in (a) and (b). 

Introducing as parameters the potential energy Vq = 0.25 cV (i.e. Vi = 0.5 
eV) and the effective electron mass m* = O.OTrrig typical for GaAs/GaAlAs 
double barrier heterostructures (see e.g. [22]), the temporal units get approx- 
imately 1.3 fs, the spatial units approximately 15 A. Therefore the width of 
the model heterostructure in figure 2 is in the range of 200 A, which is easily 
achieved by the ultrathin layers of modern semiconductor devices. 
With the aid of formulae (5) and (6) of proposition 1 the computational effort 
involved with the calculation of transmission and reflection times was reduced 
to about 5% of that corresponding to the computation by means of trajectory 
sampling. 
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Appendix A. Proof of equality (7) 



Recall that the scattering solutions under consideration are given by ^(t, •) = 
r2™0(t, •), with (/){t, •) = e^*^"*/''(/)o a solution to the free Schrodinger equation. 
The Fourier transform ip := J-'{4>o) is assumed to be localized exclusively on the 
positive half line, for which reason (j>{t, •) will be further and further localized 
to the left for large negative times (in the sense of the L^-norm || • ||). 

Figure 7 then illustrates, that for the scattering wave packet ^, formula (7) 
is a plausible conjecture. In the following a rigorous proof will be given for a 
limited class of scattering solutions. 



(t,a) ^ 




(T,,a) 

Ti->— oo 
Figure 7: Closed spaee time region. 



For the closed spacetime region in figure 7 integration of the continuity equation 
assures that 

h T b T 

j{s,a)ds+ /|*(T„OI'de- I J{s,b)ds. 







© @ 

The aim is to show that the terms ® and (2) converge towards zero in the limits 
Ti ^ —oo and & — > oo, in which case formula (7) applies. 

Proposition 2. For every scattering solution ^ with the above properties, term 
Q) vanishes in the limits Ti —>■ — oo and b —>■ oo, i.e. 

M 



lim lim 



M- 



lim / |*(-A/,Opd^ = 



(8) 



Proof: First note that the inequalities 



M 



< j |*(-M,OP d^< j l*(-A/,e)l' 



dS, 



(9) 



hold. For the given scattering solution with lim ||^(t, ■ 

i— » — oo 

cording to DoUard [18], the relation 



lim 



|*(t,Or de= lim / \Ct{Mm d^ 



0(i,-)ll 



0, ac- 



holds for every measurable set A C R and with 

1/2 



With this, inequality (9) and supp{(p} C M^ one immediately proofs (8) by 



OO OQ 

hm f\^{n,0\^d^ = lim I \CrAM(OY 

Ti~>-~OQ J Ti — > — OQ J 



d^ 



= lim 



hn 



V 



hr,. 



d^ 



lim 



|(^(fc)|^ dk 



\Lp{k)\^ dfc = 0. 



Proposition 3. For a freely evolving wave packet <f>{t, •) = e~'^^°^'^(f>o with 
Fourier transform ip := T{4>q), ip G C"'^(M) and supp{tf} C [ai, 02], < ai < 02 
(i.e. Lp 6 Cq(R-)-)J term (2) vanishes in the limits Ti — ^• —00 and b —^ 00, i.e., 



lim lim ( 



lim / j(s,M)ds^O. 

A/-too J 
-M 



(10) 



Proof: Equation (10) can immediately be shown by a stationary phase argu- 
ment. In the following the parameter reduced notation h = m = 1 will be 
used. 

The solution of the free Schrodinger equation can be written in the form 

(j>it,x) = T-^ (e-'"*^(0o)) {x) = (27r)-i/2 /" ^^ e^i^^'^im ^{j,) 



with the function cj : M ^ M, fc i-^ /c^/2 being in C°°(R). The stationary phase 
argument then states (cf, e.g., [17], appendix 1 to XI. 3) that for every open set 
A D [ai, 02] ^ {u!'{k)\k G supp{(p}}, there is a constant C > such that for all 
(i, x) e M^ with f ^ A 

|(/)(t,x)| <C(l + |x| + |f|)-\ 

The same considerations then deliver a second constant C" > such that for all 

(i, x) e r2 ^itj^ I ^ ^ 



9a; 



4>{t,x] 



dfc (i fc) ■ e'('=^-'^('='* V(fc) 



<C'(l + |a;| + |t|)-\ 



10 



We choose without loss of generality A := ] ^, 2 • 02 [. Then for a fixed t S K 
there is a M > such that for all t < t and all x > M: j ^ A (set, e.g. 
M > 2a2T). Then, Vx > M and V< < r 



\j{t,x)\ < m,x)\ 
Therefore 



d_ 
dx 



(j}{t,x) 



<c-c'-(i + |.T| + |ti: 



lim 



j{t,M) dt 



M 



< lim 

A/— foo 



\j{t,M)\dt 



-M 



<C -C hm 

A/^oc 

= C -C hm 

A/^oc 



(1 + |A/| + |i|)-2 dt 



M 



1 + sgnir) 
1 + lAf I 



1 



1 + 2 Af 1 



sgnir) 



M 



= 0. 



Now consider scattering solutions 

00 
^{t,x) = f dkip{k)4>"'{k;x)e-'''''/'^ 


with (jf^ik] x) a solution to the corresponding Lippman-Schwinger equation and 
if again exclusively localized on the positive half-line. If the ^™(fc;a;) have the 
form 



(/)"(fc;a;)=r(fc)e*'=^ 



(11) 



for a; > i? for some i? > 0, and T G C^(R) (e.g., the potential barrier), then 
the above line of reasoning applies analogously. Clearly expression (11) is only 
valid for scattering potentials with support bounded from the right. 

Appendix B. Proof of proposition 1 

By 7a: : M — > R, f i— > 7^; (t) , the integral curve of the Bohmian velocity vector field 
with initial datum x is denoted. The intervals [at,6t] := {x e M/7a;(i) S [a,^]} 
represent the initial data on configuration space, which are projected onto the 
interval [a, h\ at time t along their integral curves. Let ccc S K be the initial 
condition of the trajectory, which separates the transmitted from the reflected 
ensemble. Then the transmission time (2) gets 
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T"/ 



{tt) = da^{Q,0\' dtx[aM{ldi)) 



dtj di\^{0,O?X[a,Mm 

Ti Xc 

dt yda*(o,c)p 

[rijTf] [at.bt]n[xc.oQ[ 

|2 



dt J dei*(t,or0(e-7x.(i)). 

[Ti.Tj] [a,b] 

Analogously, the reflection time reads as 

{tb) = I dt J da^it,OfQil.At)-0■ 
[T,,Tf] [a,b] 

Now the right hand side of (5) together with (3) and (7) gets 
/ dt min<^ / \^{t,Ofd^, / \^{t,Ofd^ 

J-Ti [ [Ja "'7x0 (*) 

-minJ/ \^{t,0\^d^,f \^{t,^)\^dA 

[■Jb Jl.At) J. 

-/ max{ 



dt 

max{a, 7x^(4)} 

i 7max{a,7a;^(t)} 

2 



b.7x,(t)} 



*(i,e)i'rfc 



di / |vi'(i,e)re(e-7x.(i))rf^=<TT>. 

Analogously the right hand side of (6) together with (3) and (7) becomes 
/ dt max<^ / |vI/(i,C)pde, / \^{t,0?d^\ 

-maxi / \^{t,0\^dt / |*(t,OP'^C 






dt / I'T,^^ tM2 

' dt f \^it,o\^ ■ QhxM - d( ^< TR > 
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which completes the proof. 
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